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TECHNICAL NOTE 3120 


SPAN LOAD DISTRIBUTIONS RESULTING PROM CONSTANT VERTICAL 
ACCELERATION FOR THIN SWEPTBACK TAPERED 
WINGS WITH STHEAMWISE TIPS 


SUPERSONIC LEADING AND TR AILIN G EDGES 
By Isabella J. Cole and Kenneth Margolis 


SUMMARY 


On the basis of the linearized supersonic-flow theory, equations 
for the span load distribution resulting from constant vertical accel- 
eration (that is, linear variation of angle of attack with time) are 
derived for a series of thin sweptback tapered wings with streamwise 
tips. The analysis is valid at Mach numbers for which the wing leading 
and trailing edges are supersonic. A minor restriction is that the Mach 
line from the leading edge of either wing tip may not intersect the 
remote half wing. 

The computational results of the investigation are presented In a 
\ series of charts from which the span loadings may be obtained for given 

values of aspect ratio, taper ratio, leading-edge sweepback, and Mach 
number. For Illustrative purposes, variations of the spanwise distri- 
bution of circulation (which is proportional to the span load distribu- 
tion) with several plan-form parameters and Mach number are shown, in 
addition to some typical chordvise and spanwise pressure distributions. 


INTRODUCTION 


The aerodynamicist requires detailed information on the load dis- 
tribution over the component surfaces of an airframe. One of the most 
important considerations is the distribution of load along the wing 
span. This information can be used directly to obtain the forces, and 
moments acting on the wing itself and to estimate roughly the load on 
isolated vertical tails (for corresponding motions) and isolated hori- 
zontal tails j in addition, knowledge of the span load distribution is 
a prime requirement for the solution of problems relating to loads and 
aeroelasticity and for flow- field and other aerodynamic calculations. 
Thus, much effort has been devoted to developing methods of calculation 
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and utilizing these methods to obtain detailed load information through- 
out the range of flight speeds for wings of various plan forms undergoing 
several types of motion. 

Some of the more recent contributions to the literature on span 
load distributions at supersonic speeds are references 1 to The 
general plan form considered in these references has arbitrary aspect 
ratio, taper ratio, and sweepback; the wing tips are parallel to the 
axis of wing symmetry. At Mach numbers for which the wing leading edge 
is subsonic, equations and charts for the span load distributions 
resulting from constant angle of attack, steady rolling, steady pitching, 
and constant vertical acceleration are given in reference 1; the case of 
constant sideslip is treated in reference 2. At Mach numbers for which 
the wing leading edge is supersonic, corresponding information is given 
for constant angle of attack, steady rolling, and steady pitching in 
reference 3; calculations for the sideslip motion are given in refer- 
ence 4 . (Restrictions common to references 1 to 4 inclusive are that 
the wing trailing edge is supersonic and that the Mach line from the 
leading edge of either tip may not intersect the remote half wing . ) The 
span loading due to constant vertical acceleration for the supersonic- 
leading-edge condition has not been previously considered; the present 
paper contributes this information. 

Equations are derived herein for the spanwise distribution of circu- 
lation (which is proportional to the span load) due to constant vertical 
acceleration, that is, linear variation of angle of attack with time. 

The type of wing plan form considered is that used in references 1 to 4 ; 
namely, wings of arbitrary aspect ratio, taper ratio, and sweepback, 
and with wing tips that may be termed "streamwise" for this motion. The 
analysis is valid at supersonic speeds for which the wing leading and 
trailing edges are supersonic, provided the Mach line from the leading 
edge of one tip does not intersect the opposite half wing. 

The numerical results of the investigation are presented in a series 
of charts from which fairly rapid estimations of the load distribution 
may be obtained for given values of the parameters aspect ratio, taper 
ratio, leading-edge sweepback, and Mach number. Several illustrative 
variations of the load distribution with Mach number and wing geometry 
as well as some typical chordwise and spanwise pressure distributions 
are also presented. 


SYMBOIS 


x,y 


v 


Cartesian coordinates (see fig. l) 
free-stream velocity 
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p 

M 

M- 

B 

b 

c r 

e 

g 

h 


density of air 

free- stream Mach number 

Mach angle 

cotangent of Mach angle, y^P-1 
wing span 
root chord 

spanwise coordinate of intersection of trailing edge of wing 
and Mach line from wing apex 

spanwise coordinate of intersection of trailing edge of wing 
and Mach line reflected from wing tip 

spanwise coordinate of intersection of trailing edge of wing 
and Mach line from leading edge of wing tip 


taper ratio, 
aspect ratio, 


Tip chord 
Root chord 

2b 


c r (l + A) 
angle of sweep (see fig. l) 


m — cot At;e 


k = 


a 


AB(l + A) 


a 


AP 


COt Ajjrg 

cot Axe AB(l + a) - 4mB(l - A) 

angle of attack 

steady pitching velocity about the y-axis, positive as shown 
in figure 1 

time 

rate of change of a, with time 0^; positive a, indicates 
an acceleration downward 

difference due to constant vertical acceleration between upper- 
and lower-surface pressures; positive in sense of lift 


(AP ) a= 3_ difference due to unit angle of attack between upper- and 
lower-surface pressures; positive in sense of lift 
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(£P)q_]_ difference due to unit pitching velocity about the y-axis 
between upper- and lower- surf ace pressures; positive in 
sense of lift 

0 perturbation velocity potential due to constant vertical 

acceleration, evaluated on upper surface of wing 

0 a =l perturbation velocity potential due to unit angle of attack, 

evaluated on upper surface of wing 

0q=l perturbation velocity potential due to unit pitching velocity 

about y-axis, evaluated on upper surface of wing 

T circulation due to constant vertical acceleration at any span- 

wise station y, defined by equation ( 3 ) 

Fi, Fg components of circulation due to constant vertical accelera- 
tion; r = Ti + f 2 


circulation due to angle of attack 
r q circulation due to steady pitching velocity about y-axis 

Subscripts: 

IE refers to leading edge or evaluation along leading edge 

TE refers to trailing edge or evaluation along trailing edge 

All angles are measured in radians unless otherwise indicated. 


ANALYSIS 

Scope 


The motion considered in the present paper has been commonly termed 
constant vertical acceleration; it is a time dependent motion that 
involves a linear variation of angle of attack with time. Expressions 
for the span load distribution associated with such a motion are derived 
and calculations based on these expressions are presented for a gener- 
alized family of isolated wings (see fig. l). The wings are uncambered 
and have vanishingly small thickness; the plan form is of arbitrary 
sweepback and taper ratio with wing tips parallel to the axis of wing 
symmetry (streamwise tips). The analysis is carried out within the 
framework of the linearized supersonic-flow theory and is applicable 
at those supersonic speeds for which the wing leading and trailing 
edges are supersonic, subject to the relatively minor restriction that 
the Mach line from the leading edge of either wing tip does not inter- 
sect the remote half wing. 
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Derivation of Equations 

• The spam-rise distribution of circulation r for unyawed wings is 
related to the integrated chordwise pressure distribution (span loading) 
by the relationship 


r 



AP dx 


( 1 ) 


It should be mentioned that the span loading and the spanwise dis- 
tribution of circulation are sometimes used interchangeably. As is 
apparent from equation (l) , these two quantities are proportional to 
each other. Based on an application of the linearized supersonic-flow 
theory, the lifting pressure AP for the motion considered (that is, 
a positive d) has been shown previously (for example, eq.. ( 13 ) of 
ref. 5) to be expressible as follows: 


AP = 




( 2 ) 


where (AP)^ = 2. i- B the lifting pressure for positive unit pitching 
velocity about the y-axis and where (AP ) a=1 and (0) a =i are the 

lifting pressure and perturbation velocity potential, respectively, 
due to positive unit angle of attack. Equation (l) may then be 
rewritten: 


r = 


a. 

B 2 pV 



*£x 


(AP) a= i - 2p(0) a=1 


dx 


(3) 


It is known (for example, 


see ref. 3) that 



(AP) q=1 dx = 2pV(0 a=1 ) ra , 


oo 
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and 



(*p>o-i 


(5) 


Upon integration by parts of the middle term of equation (5) , utiliza- 
tion of equations (^) and ( 5 ), and combination of terns, equation ( 3 ) 
becomes 


T = 2a 


^ 2 (0q-l)r PE 

-D 


VB l 


2 xni \ r ^ 


( 6 ) 


According to equation (l) , the following relationships are valid: 


(Iq) 


q=i 



(AP)^! dx 


(T) 


and 


A 


( r a) a =i 



(AP) a=1 dx 


(8) 


It then follows from equations (4) and ( 5 ) that 

( r q)q=l = 2 (0q=l)rrE (9) 

and 




( 10 ) 
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Equation (6) may now be expressed in the following convenient manner 
(since and r a are directly proportional to q and a, respectively): 


r - a -la— - ^(b/ 2 ) 

B Bq(b/2) 2 B 


AB(l + + lfkmB 

2 B /2 

kmBAB(l + A) 


* g 

Val 


2 

V 



(0) a =l 


( 11 ) 


where xijrg has been replaced by its functional equivalent 


Xrjrg — 


B AB(l + A) + 14-mBkl 


jb/2 




AB(l + A)Bmk 


( 12 ) 


Formulas and charts for the quantities 


Bq(b/2) 2 


and 


L a, 


Va 


appearing in equation (ll) are presented in reference 3- Thus, only 


the integral 


TE 

X LE 


(0U 


o=l 


for the potential (0) a =l 


dx remains to be evaluated. Expressions 
for the various wing regions formed by the 


Mach line and plan- form boundaries may be obtained from reference 6 
(the value of a therein being replaced by unity). 


For convenience in calculations and presentation of results, equa- 
tion (ll) for the spanwise distribution of circulation T is subdivided 
into two components, Tq and such that T = Tq + ?£• The components 

are as follows: 


M 2 a(b/2) 2 

r AB(l + A)-^- + 4kmB 

r q b/2 

B 

Bq(b/2) 2 kmBAB(l + A) 



( 15 ) 
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0 P X TE 2(0) , 

r 2 = <xB(b/2 ) 2 / dx (l4) 

^ XlE BV(b/2) 


® r i r 2 

Formulas for the quantities — and =•, which are 

d(b/2) 2 l£ Bd(h/2) 2 

functions of the parameters AB, Bm, A , and the nondimensional span- 
wise coordinate -, are presented in table I for the combinations of 

Mach l in es n.nrt plan forms considered. The mathematical expressions 
required in the use of table I are given in the appendix. 


RESULTS AND DISCUSSION 


Equations have been derived in the preceding analysis for the span- 
wise distribution of circulation resulting from constant vertical accel- 
eration of isolated sweptback wings of arbitrary aspect ratio and taper 
ratio with streamwise tips. The analysis is valid at those supersonic 
speeds for which the wing leading and trailing edges are supersonic, 
subject to the additional minor restriction that the Mach line from the 
leading edge of one wing tip does not intersect the remote half wing. 

The pressure distributions obtained over the wing for the motion 
considered are of interest. For illustrative purposes some chordwise 
and spanwise pressure distributions, calculated from equation ( 2 ), are 
presented in figures 2 and 5, respectively. (The expressions for 
faff) 1 , (AP) Q _ 1 , and ( 0)^1 required for the calculations were 

obtained from tables V, III, and II, respectively, of reference J - in 
which both q and a therein were replaced by unity.) 

Calculations of the spanwise distribution of circulation have been 
made for values of AB from 5 to 20 for A = 0, 0.25, 0.50, 0-75, and 
1.0; for A = 1.0 calculations for AB = 2 have also been included. 

The range of Bm considered is from 1.0 to ». Results of the numerical 
calculations are presented in a series of charts (figs. 4 to 9 ) so that, 
for given value's of the parameters aspect ratio, taper ratio, leading- 
edge sweepback, and Mach number, the circulation components Tq and r 2 
(and hence r) may be readily obtained. An index to these charts is 
given in table II. 

In order to illustrate the variation of the spanwise distribution 
of circulation with different parameters , several specific examples have 
been chosen; the results are presented in figure 10 . 
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The results presented herein for constant vertical acceleration d 
may, of course, he combined with the results of reference 3 for steady 
pitching q to obtain an estimate of the span load distribution for a 
slowly (first-order frequency) oscillating wing, that is, q + & motion. 


CONCLUDING REMARKS 


On the basis of an application of linearized supersonic-flow theory, 
equations have been derived for the spanwise loading on isolated wings 
resulting from a linear angle-o'f -attack variation with time at super- 
sonic flight speeds. 

The types of wing considered have sweptback leading edges, either 
sweptback or sweptforward trailing edges, and tips that are parallel to 
the axis of wing symmetry (streamwLse tips). The analysis is applicable, 
in general, at those supersonic speeds for which the wing leading and 
trailing edges are supersonic. 

Computational results are presented in the form of generalized 
design curves which permit fairly rapid estimation of the spanwise dis- 
tribution of circulation (or span load distribution) for broad ranges 
of the parameters aspect ratio, taper ratio, leading-edge sweepback, 
and Mach number. 

The results obtained herein may be combined with corresponding 
results previously reported for the steady pitching case to obtain the 
span load distribution for a slowly (first-order frequency) oscillating 
wing. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 13, 1953* 
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APPENDIX 

LIST Of FUNCTIONS TO BE USED IN CONJUNCTION WITH TABIB I 


In order to avoid repetition of cumbersome mathematical empreaSlona appearing In the circulation 
equations of table I, a ay stem of capital -letter symbols has been used therein to denote various functions. 
Tba letters ami their functional equivalents are as follows: 



2[aB< 1 + 70-ifc + Uadi f AB(l + 70(1 -BVk)-£-+ !»Bk 

— ' _ „ , J • fabk t «fl * 70(1 - cos-1 _= '■ V*, + 

+ TO 2 ^ 2 ™ 2 - 1 L Vil -Bo|ab( 1 + A)(l - lc)^+ iMiBfcj 


*Bmi; 8 A 2 B 2 ( 1 


AB(l + 70(l + B S A.)-Z- + 4mBh 


[umBk + AB(1 + 7>)(l + kj^ljcos- 1 — l, l + fjL + 4l ™ B 1 x 

L BmpB(l + ?0(l + 4mBkj *kmB(B 2 i>2 - l) Q>/ 2 AB<1 + l)J 

lLp^i.. B W) + _^_JL + ^flfJ + fr_zL< k . i)fe . + i)l + 

f k AB(l + 70b/2 a2 b 2 (i + a) sJ H, (b/2) 2 L J 

-SSL-u + - e)JL + -Jgggf (A £ - 2)1 cos-1 

AB(l + X) b/2 A®b8(1 + X) 2 J _ h )_yL t >flk ^ 


b/2 AB(l + X)J 


(Equation continued on next page) 


a 


% 
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- ^ (k + 1 ) fTa - - i)7 + 

_(b/ 2 ) 2 


— — (bV* - 2 - k)-Z- + 

AB(l + X) b/2 


- 2 ) 

A 2 ^! + A) 2 


4-a + Aft) + _f2L , 

S -1 b/2 AB(I ± X) 

a*$i + XJL+ * 2 L 3. 
L b/2 A3(l + X)J/ 


J = ■ 


- 1 


y + 4mBk ~ j 2 
b/2 AB(l + xj 1 


y 4aBk 




4kmB 


AB(l + X) 


4kmB 

+ AB(l"+ x) 


k 2 !!® 


k^isB 3 

3m (b/2) 2 j 

[ 

Vm~R 

I' o ] 

JL + 

ImBk 

\ 

1 , 1 y 2 

1^2 + 

AB(l + X) 

y | 

aa (b/2) 2 


kmB 

b/ 2 J 


b/2/ 


Bm 


P-(l + k) 


[b/2 AB(1 

l > kkmB 
J AB(l + A) 
v 4kmB ”| 

0 + sti raj 


4-fanB ~~ | 

B(1 + X)J 


(b/2) 


2 Bn 


cosh" 


y 

4nBk 

b/2 

AB(l + x) 


kmB -2— 


b/2 


L - 


2 

i [ 

Ld* 

*kmB(B a m 2 - l) ' 

2k^B 2 m 2 - 1 



i - A( k + 1) P 2 - B 2 m 2 {k - 1)1 + 

(b/2) 2 L J 


fiftnV 
AB(l + A) 


-(-k + B*tu? 



rtAAW-S) 2 (B 2 m 2 - 1 )|aB(1 + X)^ 1 

^ +■ 4ku^| j^kmB + AB(1 + X)(l + 

A 2 B 2 (l + X) 2 

A 2 B 2 (l + X) 2 


coo-1 


-i-(l - k + 2mBk) + — ^ -r- - 2k(nB - l) 

b/2 AB(1 + A) 


-X_(l + k) + — !tS®E 

b/2 AB(l + X) 


2k' 


3 


lya 


12 + B m 2 k +■ 4mB + 


2B5m5k + 5m S B 8 ) + (fcBmk - 2B 2 m 2 k - 2B 3 m 3 k) 


4aBk(-12 + 4fflB + 5a 2 B 2 ) 


AB(l + A) 


)_ 


k(mB + 1) Y (b/2) 


-(l + Bmk) 


y IhuBIi 

b/2 AB(l + A) 


■ — k + 2mPk + 


3 


4mHk 

AB(l + A) 


+ k - kmBf 
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..( V * bV - ?] y -> ^ 

Al(l + X) b/2 A 2 ^! + X} 2 J 


Q 


k2mB^B2ma - 1 



li VwiW ® 

AB(l + 70 



a - . 1 - f- — - I y . - (k - l) ji - bV(]£ + + ^ , 0t + b't* 2 - 2)-^- + 

ittaBCB 2 !! 2 - i) L.| b 2,^ . j. I (a/a) a ai<i + x) a/a 


I6bW(bV . a) ^ - i){ab(i , x)^g t + AB(1 , X)U - k)^ 

a 2 b 2 (x+x ) 2 A 2 B®(l + TO 2 
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T/BIS I.- EXPRESSIONS FOR THE ffiAJWia DISTRIBUTION CF CIRCUIATIOH DUB TO CONSTANT VERTICAL ACCELERATION ( T - + T z ) - Concluded 

[jormilaa for the paraaetera D, S, F, O, . . . W 1 are given in appendix] 
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TABLE II.- INDEX TO CHARTS FOR SPANWTSE 


DISTRIBUTION OF CIRCULATION 


Figure 


!*• to 20 
3 
k 
6 
8 
12 
20 


1.0 to 5.0 

1.0 to 3-0 

1.0 to 00 

1.0 to 00 

1.0 to 00 

1.0 to 00 

1.0 to oo 


1.0 to 7-0 

1.0 to oo 

1.0 to 00 

1.0 to oo 

1.0 to oo 

1.0 to oo 


3 

1.0 to oo 

8(a) 

k- 

1.0 to oo 

8(b) 

6 

1.0 to oo 

8(c) 

8 

1.0 to oo 

8(d) 

12 

1.0 t6 oo 

8(e) 

20 

1.0 to oo 

8(f) 


3 

1.0 to oo 

7(a) 

k 

1.0 to oo 

7(b) 

6 

1.0 to oo 

7(c) 

8 

1.0 to oo 

7(d) 

12 

1.0 tO oo 

7(e) 

20 

1.0 to oo 

7(f) 


Figure 


2 

1.0 to oo 

9(a) 

3 

1.0 to oo 

9(b) 

k 

1.0 to oo 

9(c) 

6 

1.0 to oo 

9(d) 

8 

1.0 tO oo 

9(e) 

12 

1.0 tO oo 

9(f) 

20 

1.0 tO oo 

9(g) 
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(a) Circulation component r 1 . 



(b) Circulation component r 2 » 

Figure k.~ Distribution of circulation along span for delta wings. 

r = ii + r 2 - 
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(d) KB = 8; ~h = 0. 
Figure 5-- Continued 







(e) AB = 12; X = 0. 
Figure 5-- Continued. 





(f) AB = 20; A = 0. 
Figure 5*- Concluded. 
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b/a 

(a) AB = 3j A = 0.25. 


Figure 6.- Distribution of circulation along span for wings with 7\ = 0.25. 

r= r x + r 2 . 
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b/z 

(f) AB = 20; ~h = 0.25. 
Figure 6.- Continued. 
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(To) AB = k} A = 0 . 50 . 
Figure 7 .- Continued. 





(d) AB = 8; * = 0.50. 
Figure 7 *- Continued. 
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bjt 

(e) AB = 12; ~K = 0.50. 


Figure 7*- Continued. 


(e) AB = 12 ; = 0.50. Concluded 

Figure 7 .- Continued. 



41 




~&J£ 


(f) AB = 20; X = 0.50. 
Figure 7 .- Coni 


inued . 


(f) AB = 20; Tv = 0.50. Concluded 
Figure 7-" Concluded. 
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(c) AB = 6; \ - 0. 75- 

Figure 8.- Continued. 











(e) AB = 12; A = 0.75. 
Figure 8.- Continued. 
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£>/a 


(b) AB = 3; A = 1.0. 
Figure 9 -~ Continued. 
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(c) A3 = k; A = 1.0. 
Figure 9-- Continued. 
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HJ2. 

(g) AB = 20; A = 1.0. 


Figure 9*- Continued. 






(g) AB = 20; 7 \ = 1.0. Concluded. 
Figure 9.- Concluded. 
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(a) Variation with Mach number. A = k; A = 30 


(b) Variation with aspect ratio. M = 1.53; A = 50 ; A = O.pO. 

Figure 10.- Some examples illustrating the variations of distribution of 
circulation along span with Mach number, aspect ratio, sweepback, and 
taper ratio. 










